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Recent advances in synthesis and electrical characterization of nanofibers and nanotubes made out 
of various conjugated polymers attract attention of the research community to studies of transport 
properties of these materials. In this work we present a theoretical analysis of electron transport in 
polymer nanofibers assuming them to be in conducting state. We treat a conducting polymer as a 
network of metallic-like grains embedded in poorly conducting environment, which consists of ran- 
domly distributed polymeric chains. We analyze the contribution from intergrain electron resonance 
tunneling via intermediate states localized on the polymeric chains between the grains. Correspond- 
ingly, we apply the quantum theory of conduction in mesoscopic systems to analyze this transport 
■ mechanism. We show that the contribution of resonance electron tunneling to the intergrain elec- 

tron transport may be predominating, as follows from experiments on the electrical characterization 
of single polyaniline nanofibers. We study the effect of temperature on the transport characteris- 
tics. We represent the thermal environment as a phonon bath coupled to the intermediate state, 
which provides electron tunneling between the metallic-like grains. Using the Buttiker model within 
the scattering matrix formalism combined with the nonequilibrium Green's functions technique, we 
show that temperature dependencies of both current and conductance associated with the intergrain 
electron tunneling, differ from those typical for other conduction mechanisms in conducting poly- 
mers. Also, we demonstrate that under certain conditions the phonon bath may cause suppression 
of the original intermediate state accompanied by emergence of new states for electron tunneling. 
The temperature dependencies of the magnitudes of the peaks in the transmission corresponding to 
"j^ ' these new states are analyzed. 
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During the past decade, transport properties of conducting polymers such as doped polyacetylene and polyaniline- 
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polyethylene oxides, were and still are intensively studied [l|, |2|. These materials are significant mostly due to 
various possible applications in fabrication of nanodevices. Polymer-based devices should have advantages of low 
cost and flexible, controlled chemistry. Also, there are some unresolved problems concerning the physical nature of 
, charge transfer mechanisms in conducting polymers, which make them interesting subjects for fundamental research. 
' Chemically doped polymers are known to be very inhomogeneous. In some regions polymer chains are disorderly 
I/"") , arranged, forming an amorphous, poorly conducting substance. In other places the chains are ordered and densely 
packed [I i]. These regions could behave as metallic-like grains embedded in the disordered environment. The 
t-H , fraction of metallic-like islands in bulk polymers varies depending on the details of the synthesis process. In practical 
samples such islands always remain separated by disordered regions and do not have direct contacts. In some cases, 
electronic states are delocalized over the grains, and electrons behave as conduction electrons in conventional metals. 
In these cases electrons motion inside the grains is diffusive with the diffusion coefficient D = \v 2 f t (vf is the Fermi 
velocity, and r is the scattering time). In whole, electron transport in conducting polymers shows both metallic and 
nonmetallic features, and various transport mechanisms contribute to the resulting pattern. 

An important contribution to the conduction in these substances is provided by the phonon-assisted electron 
hopping between the conducting islands and/or variable range hopping between localized electronic states. The effect 
of these transport mechanisms strongly depends on the intensity of stochastic nuclear motions. The latter increases as 
temperature rises, and this brings a significant enhancement of the corresponding contributions to the conductivity. 
The temperature dependence of the "hopping" conductivity a(T) is given by the Mott's expression [Bjj: 

( 7(T)=a(0)exp[-(T /Tf] (1) 

where Tq is the characteristic temperature of a particular material, and the parameter p takes on values 0.25, 0.33 
or 0.5 depending on the dimensions of the hopping processes. Also, it was suggested that phonon-assisted transport 
in low-dimensional structures such as nanofibers and nanotubes, may be substantially influenced due to electron 
interactions [1,01 • This results in the power-low temperature dependencies of the conductance G(T) at low values 
of the bias voltage V (eV < kT, k being the Boltzmann constant), namely: G ~ T a . Experimental data for the 
conductance of some nanofibers and nanotubes match this power-low reasonably well, bearing in mind that the value 
of the exponent a varies within a broad range. For instance, a was reported to accept values about 0.35 for 
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carbon nanotubes and a ~ 2.2 -f- 7.2 for various polyacetylene nanofibers [§-11 . In general, hopping transport 
is very important in disordered materials with localized states. For this kind of transport phonons play part of 
a source of electrical conductivity. Accordingly, the hopping contribution to the conductivity always increases as 
temperature rises, and more available phonons appear. When polymers are in the insulating state, the hopping 
transport predominates and determines the temperature dependencies of transport characteristics. 

In conducting state of conducting polymers free charge carriers appear, and their motion strongly contributes to 
the conductance. While moving, the charge carriers undergo scattering by phonons and impurities. This results in 
the conductivity stepping down. Metallic-like features in the temperature dependencies of dc conductivity of some 
polymeric materials and carbon nanotubes were repeatedly reported. For instance, the decrease in the conductivity 
upon heating was observed in polyaniline nanofibers in Refs. fl2| and [l3j ]. However, this electron diffusion is 
not a unique transport mechanism responsible for the occurrence of metallic-like behavior in the dc conductivity 
of conducting polymers. Prigodin and Epstein suggested that the electron tunneling between the grains through 
intermediate resonance states on the polymer chains connecting them, strongly contributes to the electron transport 
[Hj]. This approach was employed to build up a theory of electron transport in polyaniline based nanofibers [HI 
providing good agreement with the previous transport experiments fl6j | . Considering the electron tunneling through 
the intermediate state as a mechanism for the intergrain transport, we see a similarity between the latter and electron 
transport mechanisms typical for tunnel molecular junctions. In the case of polymers, metallic-like domains take on 
part of the leads, and the molecular bridge in between is simulated by intermediate sites. The effect of phonons on 
this kind of electron transport may be very significant. These phonons bring an inelastic component to the intergrain 
current and underlie the interplay between the elastic transport by the electron tunneling and the thermally assisted 
dissipative transport. Also, they may cause some other effects, as was shown while developing the theory of conduction 
through molecules 



II. Electron-resonance tunneling as a transport mechanism in conducting polymers 

Here, we concentrate on the analysis of the electric current-voltage characteristics and conductance associated 
with the resonance tunneling transport mechanism. Considering the electron intergrain resonance tunneling, the 
transmission coefficient is determined with the probability of finding the resonance state in between the grains. 
The latter is estimated as P ~ exp(— £/£) (L is the average distance between the adjacent grains, and £ is the 
localization length for electrons), and it takes values much greater than the transmission probability for sequental 
hoppings along the chains, Ph ~ exp(— 2L/£) [3]. Nevertheless, the probability for existence of a resonance state 
at a certain chain is rather low, so only a few out of the whole set of chains connecting two grains are participating 
in the intergrain electron transport. Therefore, one could assume that any two metallic domains are connected by a 
single chain providing an intermediate state for the resonance tunneling. All remaining chains can be neglected for 
they poorly contribute to the transport compared to the resonance chain. Within this approximation the "bridge" 
linking two islands is reduced to a single electron state. Realistic polymer nanofibers have diameters within the range 
20 -7- 100 nm, and lengths of the order of a few microns. This is much greater than a typical size of both metallic-like 



grains and intergrain separations, which take on values ^ 5 ™ 10 nm (see e.g. Refs. 16| and [25|). Therefore, we may 
treat a nanofiber as a set of working channels connected in parallel, any single channel being a sequence of grains 
connected with the resonance polymer chains. The net current in the fiber is the sum of currents flowing in these 
channels, and the voltage V applied across the whole fiber is distributed among sequental pairs of grains along a 
single channel. So, the voltage AV applied across two adjacent grains could be roughly estimated as AV ~ VL/Lq 
where L is the average separation between the grains, and Lq is the fiber length. In practical fibers the ratio AV/V 
may take on values of the order of 10~ 2 -j- 10~ 3 . 

In further current calculations we treat the grains as free electron reservoirs in thermal equilibrium. This assumption 
is justified when the intermediate state (the bridge) is weakly coupled to the leads, and conduction is much smaller 
than the quantum conductance Go = 2e 2 /h (e, h are the electron charge, and the Planck constant, respectively). Due 
to the low probabilities for the resonance tunneling between the metallic islands in conducting polymers, the above 
assumption may be considered as a reasonable one. So, we can employ the well-known expression for the electron 
current through the molecular junction (2rj |. and we write: 

I=™ dET{E)[h{E)-f 2 {E)]. (2) 



Here, n is the number of the working channels in the fiber, fi^(E) are Fermi functions taken with the different 
contact chemical potentials for the grains. The chemical potentials differ due to the bias voltage AV applied 
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FIG. 1: Schematic drawing illustrating the intergrain electron transport in the presence of dissipation [291 ]. Rectangles 
correspond to the barriers separating the adjacent metallic-like islands from the intermediate state (the bridge), and the 
triange stands for a scatterer attached to the bridge. 



across the grains: 

Hx = E F + (1 - rj)eAV; /jl 2 = E F - TjeAV. (3) 

In these expressions ([2]), (j3|) the parameter r/ characterizes how the voltage AV is divided between the grains, Ep 
is the equilibrium Fermi energy of the system including the pair of grains and the resonance chain in between, and 
T(E) is the electron transmission function. 

The general approach to the electron transport studies in the presence of dissipation is the reduced dynamics density- 



matrix formalism (see, e.g., Refs. [27| and |28[). This microscopic computational approach has the advantages of 
being capable of providing the detailed dynamics information. However, this information is usually more redundant 
than necessary, as far as standard transport experiments in conducting polymer nanofibers are concerned. There 
exists an alternative approach using the scattering matrix formalism and the phenomenological Buttiker dephasing 
model [29[ . Adopting this phenomenological model we are able to analytically study the problem, and the results 
agree with those obtained by means of more sophisticated computational methods, as was demonstrated in the earlier 
works [3fjj . 

Within the Buttiker model we treat the intergrain electron transport as a multichannel scattering problem. In the 
considered case the "bridge" between two adjacent grains inserts a single electron state. Therefore, an electron could 
be injected into the system (including two metallic- like domains and the intermediate "bridge" in between) and/or 
leave from there via four channels presented in the Fig. 1. The electron transport is a combination of tunneling 
through two barriers (the first one separates the left metallic island from the intermediate state and the second 
separates this state from the right island, supposing the transport from the left to the right). Inelastic effects are 
accounted for by means of a dissipative electron reservoir attached to the bridge site. The dissipation strength is 
characterized by a phenomenological parameter e, which can take values within the range [0, 1]. When e = the 
reservoir is detached from the bridge, which corresponds to the elastic and coherent electron transport. The greater 
is e value the stronger is the dissipation. In the Fig. 1 the barriers are represented by the squares, and the triangle 
in between imitates a scatterer coupling the bridge to a dissipative electron reservoir. 

Incomin g pa rticle fluxes (Jj) are related to those outgoing from the system (Jj) by means of the transmission 
matrix T [23, 0: 



•/; }2'j->- i<m< 4 - ( 4 ) 

i 

Off-diagonal matrix elements Tji(E) are probabilities for an electron to be transmitted from the channel i to the 
channel j, whereas diagonal matrix elements Tu(E) are probabilities for its reflection back to the channel i. To 
provide charge conservation, the net particle flux in the channels connecting the system with the reservoir must be 
zero. So we have: 

j 3 + Ji-4-4 = o. (5) 

The transmission function T(E) relates the particle flux outgoing from the channel 2 to the flux incoming to the 
channel 1, namely: 



f 2 = T{E)J x . 



(6) 
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Using Eqs. (0| and ([5]) we can express the transmission function in terms of the matrix elements T^. The latter are 
related to matrix elements of the scattering matrix S, which expresses the outgoing wave amplitudes b[ , b' 2 ■ , as 



linear combinations of the incident ones b\, & 2 , 03, a a ■ Tij — \Sij\ 2 . In the considered case of a single site bridge the 



S matrix takes the form lf| 30| : 
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where Z = 1 — a 2 rir 2 , a = y/1 — e, /3 = ^/e, ri )2 and ii. 2 are the transmission an reflection coefficients for the 
barriers (|ti, 2 | 2 + \r ia \ 2 = 1). 

When the bridge is detached from the dissipative reservoir T(E) = |«Si 2 | 2 . On the other hand, in this case we can 
employ a simple analytical expression for the electron transmission function [311 ]: 

T(E)^4A 1 (E)A 2 (E)\G(E)\ 2 1 (8) 

where Ai j2 (i?) = — ImSi i2 (£'). In this expression, self-energy terms Si j2 appear due to the coupling of the metallic- 
like grains to the intermediate state (the bridge). The retarded Green's function for a single-site bridge could be 
approximated as follows: 

G ^ - E~h^ (9) 

where E\ is the site energy. The width of the resonance level between the grains is determined by the parameter 
T = Ai + A 2 + T en (T en describes the effect of energy dissipation). Further we consider dissipative effects originating 
from electron-phonon interactions, so, T en is identified with T p h- 

Equating the expression (jHJl and |Si 2 | 2 we arrive at the following expressions for the tunneling parameters 5\^{E) : 



S 12 (E)^t 2 12 = - 2Al2 (10) 
1,2 ^{E-E^+T 2 



Using this result we easily derive the general expression for the electron transmission function: 

_ g(E)(l + a 2 )[g{E){\ + a 2 ) + l~ a 2 ] 
1 ' [g(E)(l-a 2 ) + l + a 2 ] 2 [ ' 

where: 



g(E)=2j 1 o- (12) 

■ n ' y (E- e x ) 2 + r 2 v ; 

To simplify further analysis we approximate the self-energy terms Ai j2 as constants: Ai, 2 « W 2 2 /<Ji } 2- Here, 
W\ t 2 are coupling strengths characterizing the coupling of the grains to the bridge and 7i j2 characterize interatomic 
couplings inside the grains (leads). Simulating the leads by semiinfinite chains of identical sites, as was first suggested 
by D'Amato and Pastawski [32], one may treat the parameters 7^2 as coupling strengths between the nearest sites 
in these chains. Then the elastic electron transmission (e = 0) shows a sharp peak at the energy E\ (see Fig. 2a), 
which gives rise to a steplike form of the volt-ampere curve presented in Fig. 2b. When the reservoir is attached to 
the electronic bridge (e ^ 0) , the peak in the transmission is eroded. The greater is the value of the parameter e , the 
stronger is the erosion. When e takes on the value of 0.7 the peak in the electron transmission function is completely 
washed out as well as the steplike shape of the I — V curve. The latter becomes linear, corroborating the well-known 
Ohmic law for the sequential hopping mechanism. So, we see that the electron transmission is affected by stochastic 
nuclear motions in the environment of the resonance state. When the dissipation is strong (e.g. within the strong 
thermal coupling limit), the inelastic (hopping) contribution to the intergrain current predominates, replacing the 
coherent clastic tunneling. Typically, at room temperatures the intergrain electron transport in conducting polymers 
occurs within an intermediate regime, when both elastic and inelastic contributions to the electron transmission are 
manifested. 

The described approach was succesfully employed to analyze experimental results on the electrical characterization 
of polyaniline-polyethylene oxides (PANi) nanofibers reported by Zhou et al [l6j]j and Pinto et al !33J. In these 
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FIG. 2: Dissipation effect on the effective transmission function (a) and the current voltage characteristics (b). The curves are 
plotted at T = 70K, Ai = A2 = WmeV, e = 0.0, 0.1, 0.3, 0.7 from the top to the bottom. Calculated (dash-dot lines) and 
experimental (solid lines) current (nA) - voltage (V) characteristics for PANi nanofibers: (c) samples 1 (SI), and 2 (S2), and 
(d) sample 3 (S3). 



experiments single-fiber electrical characterization was carried out at T ~ 300K. The experiments revealed that 
current-voltage curves for conducting nanofibers were non-Ohmic. Conducting samples included a 70-nm-diameter 
nanofiber (sample 1), a pair of 18-nm- and 25-nm-diameter fibers connected in parallel (sample 2) and a fiber whose 
diameter was steeply reduced in the middle from 70 nm to 20 nm (sample 3). The latter did exhibit a very asymmetric 
rectifying current-volltage characteristic. The expression ([2]) was employed to compute the current originating from 
electron resonance tunneling in the nanofibers |15| . To evaluate the number of working channels in a nanofiber it 
was assumed that for certain doping and crystallinity rates the number of channels is proportional to the fiber cross- 
sectional area. Also, it was taken into account that the contact of the fiber surfaces with atmospheric gases reduced 
the number of working channels. Accepting the value 5 nm for both average grain size and intergrain distance we 
estimated that the 70 nm fiber could include about 30 — 40 conducting channels, and we could not expect more than 
two working channels in the pair of fibers (sample 2). This gives the ratio of conductions ~ 15 — 20, which is greater 
than the ratio of cross-sectional areas of the samples A\ jA% w 5 . The difference originates from the stronger dedoping 
of thinner fibers. Calculating the current, appropriate values for Wx t 2 providing the best match between calculated 
and experimental I — V curves were found applying the least-squares procedure. We estimated values of the coupling 
strengths W\ t 2 in PANi fibers as 2.0 — 3.5meV. As shown in the Fig. 2, computed volt-ampere characteristics match 
the curves obtained in the experiments reasonably well. The agreement between the presented theoretical results and 
experimental evidence proves that electron tunneling between metallic-like grains through intermediate states at the 
resonance chains can really play a significant part in the transport in conducting polymers provided that they are in 
the metallic state. 

However, we remark that the good fitting between the theory and the experiment demonstrated in the work 15 1 
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was achieved assuming that dissipation was very weak (a « 0.05). This assumption could hardly be justified while 
one is observing electron transport in nanofibers at room temperature. So, the phenomenological approach employed 
in this section has some significant shortcomings. Its main disadvantage is that the dissipative effects are described 
in terms of a phenomenological parameter e, whose dependence of characteristic factors affecting the transport (such 
as temperature, electron-phonon coupling strength and some others) remains unclear. It is necessary to modify the 
Buttiker's model to elucidate the relation of the parameter e to the relevant energies characterizing electron transport 
in the considered systems. To obtain the desired expression for e , we compare our results with those presented in Ref. 



17j . In that work the inelastic correction SI to the coherent tunnel current via a single-site bridge is calculated using 
the nonequilibrium Green's functions formalism. The relevant result is derived in the limit of weak electron-phonon 
interaction when r p ^ <§C A 12 . It is natural to assume that the dissipation strength e is small within this limit, so 
we expand our expression for T(E) in powers of e. Keeping two first terms in this expansion, and assuming that 
Ai = A 2 = A we obtain: 

T{E)^g\E)(l+ t -[l-2g\E)]). (13) 

We employ this approximation to calculate the current through the bridge, and we arrive at the following expression 
for 81: 

e f°° eA {E- Erf -4A 2 



/ A 
dE-^—p el (E) 
-00 1 +Lph 



(£-£i) 2 +4A 2 

x [h(E)-f 2 (E)}. (14) 
Here, p e i(E) is the electron density of states at the bridge: 

Pel (E) = --lmG(E). (15) 

TT 



Comparing the expression (|14[) with the corresponding result of Ref. 17j, we find that these two are consistent, 



and we get [34j |: 



When A 3> T p h (e < 1) the bridge coupling to the dissipative reservoir is weak, and the elastic electron tunneling 
predominates. The opposite limit A <C T p h (e ~ 1) corresponds to the completely incoherent phonon-assistcd electron 
transport. 

III. Temperature Dependencies of electron transport characteristics in conducting polymer 

nanofibers 

Now, we turn to studies of temperature dependencies of the electron current and conductance associated with the 
intcrgrain electron tunneling. It is known that various conduction mechanisms may simultaneously contribute to 
the charge transport in conducting polymers, and their relative effects could significantly differ depending on the 
specifics of synthesis and processing of polymeric materials. The temperature dependencies of the resulting transport 
characteristics may help to identify the predominating transport mechanism for a paricular sample under particular 
conditions. The issue is of a significant importance because the relevant transport experiments are often implemented 
at room temperature, so that the influence of phonons cannot be disregarded. Therefore, we study the effect of 
temperature (stochastic nuclear motions) on the resonance electron tunneling between metallic-like grains (islands) 
in polymer nanofibers. In these studies we consider the dissipative reservoir attached to the "bridge" site as a phonon 
bath, and we assume that the phonon bath is characterized by the continuous spectral density J(w) of the form [35| : 

J(u) = Jo— exp(-— ) (17) 

where Jo describes the electron-phonon coupling strength, and lo c is the cut-off frequency of the bath, which 
determines the thermal relaxation rate of the latter. The expression for T p h(E) was derived in earlier works [17ll36l .l37| 
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using nonequilibrium Green's functions approach. Using these results and the expression (1171) , we may present T p h (E) 
as follows: 

f uj ( UJ 

T p h =2ttJq duj — exp 

J u c \ UJ C 

x {N(uj)[p el (E- hu) + Pe i{E + huj)} 

+ [1 - n(E - haj)]p e i(E - hw) 

+ n{E + tmj)p el (E + hu))}. (18) 

Here, 

and N(uj) is the Bose-Einstein distribution function for the phonons at the temperature T. The asymptotic expression 
for the self-energy term T p h depends on the relation between two characteristic energies, namely: huj c and kT (k is 
the Boltzmann constant). At moderately low or room temperatures kT ~ 10 -j- 30meV. This is significantly greater 
than typical values of hoj c (huj c ~ IraeV fl7|). Therefore, in further calculations we assume huj c <C fcT. Under 
this assumption, the main contribution to the integral over uj in the Eq. ()18|) originates from the region where 
uj <C w c <C kT/h, and we can use the following approximation [38j: 

2rA(j ,^,r) 

- (^Tp- C 2 °) 



Here, T = Ai + A 2 + T 



A = ^^l 2 ^ + 1 l (21) 



where £(2; kT/huj c + 1) is the Riemann £ function: 



C = (2;fcT/^ + l) = E (22) 

n— 1 v ' ' 



Under huj c <C kT, we may apply the estimation A AkTJo- 

Solving the equation (1201) we obtain a reasonable asymptotic expression for T p h 



T ph = ^> . ( 23 ) 

2 4(^) 2 + (i + yTT7) 2 

where p 2 = 8A/(Ai + A2) 2 . Substituting this expression into Eq. (fTB)) we arrive at the result for the dissipation 
strength e : 



e=i ^±VHZ} (24) 

2 4(^) 2 + i(i + v ^T7) 3 

This expression shows how the e depends on the temperature T, the electron-phonon coupling strength Jo, and the 
energy E. In particular, it follows from the Eq. (|24j) that e reachs its maximum at E — Ex, and the peak value of 
this parameter is given by: 



+ ■ 



The maximum value of the dissipative strength is determined with two parameters, namely, T and Jo- As illus- 
trated in the Fig. 3, e max increases when the temperature rises, and it takes on greater values when the electron- 
phonon interaction is getting stronger. This result has a clear physical sense. Also, as follows from the Eq. (|2~3| . the 
dissipation parameter exhibit a peak at E = E\ whose shape is determined by the product kT 'Jq. When either Jq 
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FIG. 3: Temperature dependencies of the maximum dissipation parameter e (left panel) and the electron transmission (right 
panel). The curves are plotted at To = 50K, Ai = A2 = imeV, E = Ex = 0, Jo = 9.0meV (solid line), Q.OmeV (dashed 
line), 3meV (dash-dotted line), l.bmeV (dotted line). 



or T or both enhance, the peak becomes higher and its width increases. The manifested energy dependence of the 
dissipation strength allows us to resolve the above mentioned difficulty occurring when the inelastic contribution to 
the electron transmission function is estimated using the simplified approximation of the parameter e as a constant. 
When the energy dependence of e is accounted for, the peak in the electron transmission at E — E\ may be still 
distinguishable when e max takes on values as big as 0.5. 

The obtained result enables us to analyze the temperature dependencies of the electric current and conductance 
of the doped polymer fibers assuming that the resonance tunneling predominates in the intergrain electron transport 
in the absence of phonons. Current-voltage characteristics and voltage dependencies of the conductance G = dl/dV 
computed using the expressions ©, (fTTj) . ([Ml) are presented in the Fig. 4. We see that at low values of the applied 
voltage the electron-phonon coupling brings an enhancement in both current and conductance, as shown in the top 
panels of the Fig. 4. The effect becomes reversed as the voltage grows above a certain value (see Fig. 4, the bottom 
panels) . This happens because the phonon induced broadening of the intermediate energy level (the bridge) assists 
the electron transport at small bias voltage. As the voltage rises, this effect is surpassed by the scattering effect of 
phonons which resists the electron transport. When the electron-phonon coupling strengthens, the TV curves lose their 
specific shape typical for the elastic tunneling through the intermediate state. They become closer to straight lines 
corresponding to the Ohmic law. At the same time the maximum in the conductance originating from the intergrain 
tunneling gets eroded. These are the obvious results discussed in some earlier works (see e.g. Ref. 30]). The relative 
strength of the electron-phonon interaction is determined by the ratio of the electron-phonon coupling constant Jo 
and the self-energy terms describing the coupling of the intermediate state (bridge) to the leads Ai 2. The effect of 
phonons on the electron transport becomes significant when Jo > Ai^. Otherwise, the coherent tunneling between 
the metallic-like islands prevails in the intergrain electron transport, and the influence of thermal phonon bath is 
weak. Again, we may remark that Jo and T are combined as fcTJo in the expression (|24p . Therefore, an increase 
in temperature at a fixed electron-phonon coupling strength enhances the inelastic contribution to the current in the 
same way as the previously discussed increase in the electron-phonon coupling. 

Now, we consider temperature dependencies of the electric current and conductance resulting from the intergrain 
electron tunneling via the intermediate localized state. These dependencies are shown in the Fig. 5. The curves in 
the figure are plotted at low bias voltage (V = 0.3V, AV/V = 0.005) and T = 50JC, so eAV < kT. This regime is 
chosen to compare the obtained temperature dependencies with those typical for the phonon assisted hopping transport 
discussed. We see that the temperature dependence of the tunnel current shown in the left panel of Fig. 5 crucially 
disagrees with the Mott's expression ([T]). The tunnel current decreases as temperatute rises being proportional to 
(To/T)' 3 , and the exponent takes on values close to unity. 

Already it was mentioned that the drop in the conductivity upon heating a sample was observed in polymers 
and carbon nanotubes. However, such metallic-like behavior could originate from various dc transport mechanisms. 
Correspondingly, the specific features of temperature dependencies of the conductivity and/or current vary depending 
on the responsible conduction mechanism [39(. The particular temperature dependence of the electron tunneling 
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FIG. 4: Current (left panels) and conductance (right panels) versus voltage. The curves are plotted at To = 50K, T = 
30K, E x = 40meV, Ai = A 2 = 4meV, n = 1, AV7V = 0.005, assuming J = 9meV (solid line), J = 6meV (dashed line), 
Jo = 3meV (dash-dotted line), and J = 1.5meV (dotted line). Go = 2e 2 /h is the quantum conductance. 



current shown in Fig. 5 differs from those occurring due to other transport mechanisms. Such dependence was observed 
in the experiment on the electron transport in a single low-defect-content carbon nanotube rope, whose metallic-like 
conductivity was manifested within a wide temperature range (T ~ 35 -j- 3Q0K), as reported by Fisher et al [io| . 
The conductivity temperature dependence observed in this work could be approximated as a (T)/o (300) ~ a + bTo/T 
where A, b are dimensionlcss constants. The approximation includes the temperature independent term, which 
corresponds to the Drude conductivity. The second term is inversely proportional to the temperature in agreement 
with the results for the current shown in the Fig. 5 (left panel). It is also likely that a similar approximation may be 
adopted to describe experimental data obtained for chlorate-doped polyacetylene samples at the temperatures below 
100K [41]. In both cases we may attribute the contribution proportional to 1/T to the resonance electron tunneling 
transport between metallic-like islands. 

The conductance due to the electron intergrain tunneling reduces when the temperature increases, as shown in the 
right panel of the Fig. 5. Irrespective of the electron-phonon coupling strength we may approximate the conductance 
by a power law G ~ T a where a takes on values close to —1. This agrees with the results for the current. At 
higher bias voltage the temperature dependence of the current changes, as shown in the Fig. 6. The curves shown 
in this panel could be approximated as log(//7o) ~ c 4- dTo/T; c,d being dimensionless constants. This resembles 
typical temperature dependencies of the tunneling current in quasi-one-dimensional metals which were predicted for 
conducting polymers being in a metal state (see e.g. Ref. 0]). 
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T n /T log(T/T n ) 



FIG. 5: Temperature dependence of current (left panel) and conductance (right panel) at low voltage bias (V = 0.3V). The 
curves are plotted at To = 50K, Jo — 9meV (solid line), Jo = 6meV (dashed line), Jo = 3meV (dash-dotted line), and 
Jo = 1.5meV (dotted line). The values of remaining parameters are the same as in the Fig. 4. The current To is computed 
using Jo = 3meV; T /T = 1.5. 



IV. Effect of phonon induced electron states on the transport properties of conducting polymer fibers 

Studies of dissipative effects in the intergrain electron transport in conducting polymers may not be restricted 
with the plain assumption of direct coupling of the bridge site to the phonon bath. Other scenarios can occur. In 
particular, analyzing electron transport in polymers, as well as in molecules, one must keep in mind that besides the 
bridge sites there always exist other nearby sites with close energies. In some cases the presence of such sites may 
strongly influence the effects of stochastic nuclear motions on the characteristics of electron transport. This may 
happen when the nearby sites somehow "screen" the bridge sites from direct interactions with the phonon bath. Here, 
we elucidate some effects which could appear in the electron transport in conducting polymer fibers in the case of 
such indirect coupling of the bridge state to the phonon bath. 

We mimic the effects of the environment by assuming that the side chain is attached to the bridge, and this chain 
is affected by phonons. This model resembles those used to analyze electron transport through macromolecules [l9j . 
The side chain is introduced to screen the resonance state making it more stable against the effect of phonons. We 
assume that electrons cannot hop along the side chain, so it may be reduced to a single site attached to the resonance 
site (the bridge). 

Within the adopted model the retarded Green's function for the bridge acquires the form (HEH: 



G^ 1 (E) = E — Ei — Ai — A 2 — w 2 P{E). (26) 

The first four terms in this expression represent the inversed Green's function for the resonance site coupled to the 
two grains, and the factor w is the hopping integral between the bridge and the attached side chain. The term P{E) 
represents the effect of the phonons and has the form: 

/>oo 

P(E) =-i I dt exp [it(E - E x + 5 + i0 + )] 
Jo 

x{(l-/)exp[-F(t)]+/exp[-F(-t)]} (27) 

with E\ being the on-site energy for the side site, which is close to the bridge site energy Ei, exp[—F(t)] being a 
dynamic bath correlation function, and / taking on values 1 and when the attached site is occupied and empty, 
respectively. 

Characterizing the phonon bath with a continuous spectral density J(w) given by Eq. (|17p one may write out the 
following expressions for the functions F(t) and 5 : 



F( t) = r p { .) 

Jo w 



1 _ e -i»t , 2[l-cos(wt)] 



exp(Hu) / kT) — 1 



(28) 
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S = / —J(u) = Jo- 



Within the short time scale (u> c t < 1) the function F(t) could be presented in the form [4J 



Jo 



hu> c I 2 



= — — { - ln[l + (uj c t) 2 } + i arctan(w c i) + K(t) 



where 



K{t) = {kTyt 2 C 2 



kT 



1 



(29) 



(30) 



(31) 



Here, C, (2; kT/hu c + l) is the Riemann £ function. The asymptotic expression for K(t) depends on the relation 
between two parameters, namely, the temperature T and the cut-off frequency lo c of the phonon bath. Assuming 

kT 3> hLj c 



In the opposite limit fkj c 3> kT we obtain: 



K(t) 



K{t) 



kT 

YlljJr 



UJ c tf 



ir 2 /kTty 



(32) 



(33) 



Also, we may roughly estimate K(t) within the intermediate range. Taking kT w Huj c we arrive at the approximation 



K(t) 



2 (kTt/h) 2 where a 2 is a dimensionless constant of the order of unity. Correspondingly, within the short 



time scale we can omit the first term in the expression (|30[) . and we get: 

TUjJr 



(34) 



where K{t) is given by either Eq. (|32p or Eq. (|3^|) depending on the relation between haj c and kT. 

Within the long time scale uj c i 3> 1, and provided that temperatures are not very low (kT 3> fiw c ), we may 
present the function K(t) as: 



K(t) 



2kTt 



dz 



(1 



cos z e 



-z/u) c t 



irkTt 



(35) 



Now, the term K(t) is the greatest addend in the expression for Fit), so the latter could be approximated as: 
F(t) ~ irkTtJo/uich 2 . The same approximation holds within the low temperature limit when uj c 3> kT /h ^> 

Using the asymptotic expression ([34|. we may calculate the contribution to P{E) coming from the short time scale 
(uj c t < 1). It has the form: 



P 1 (E) 




(36) 



where $(z) is the probability integral. When both Huj c and kT have the same order of magnitude the expression 
for P(E) still holds the form (l3l)|) . At kT <C hui c , the temperature kT in the expression (|3"6"|) is to be replaced by 
fvjj c . We remark that under the assumption kT 3> huj c the function P\{E) does not depend on the cut-off frequency 
u) c , whereas at hw c ^> kT it does not depend on temperature. The long time (ui c t ^> 1) contribution to P(E) could 
be similarly estimated as follows: 



P2(E) = 



1 



E — Ei + J + inJ kT/huj c 



(37) 



Comparing these expressions (|3"6")l and (|3"T)) we see that the ratio of the peak values of p2(E) and Pi(E) is of 
the order of {h 2 ui 2 / J^kT) 1 / 2 . Therefore, the term P\(E) predominates over p2(E) when the temperatures are 
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moderately high (hu> c < kT) and the electron-phonon interaction is not too weak Jo/Huj c ~ 1. Usually, experiments 
on the electrical characterization of conducting polymer nanofibers are carried out at T ~ 100 4- 300-ftT, so in further 
analysis we assume that (h 2 oj^/ JokT) 1 / 2 -C 1, and the term Pa(-S) could be omitted. As shown in the Fig. 6, the 
imaginary part of P{E) exhibits a dip around E — E± and the width of the latter is determined by the product of 
the temperature kT (or JUu c ) and the constant Jo characterizing the strength of the electron-phonon interaction. 
When either factor increases, the dip becomes broader and its magnitude reduces. 

The presence of the term w 2 P(E) gives rise to very significant changes in the behavior of the Green's function given 
by the Eq. (|26|) . Using the flat band approximation for the self-energy corrections and disregarding for a while all 
imaginary terms in the Eq. (j2"6| , we find that two extra poles of the Green's function emerge. Assuming Ei = E\ = 
and kT ^> hw c , these poles are situated at: 



E = ±2y/j Q kT\\n{2J kT/w 2 )\. (38) 

The poles correspond to extra electron states which appear due to electrons coupling to the thermal phonons. These 
new states are revealed in the structure of the electron transmission T(E) given by Eq. ©. The structure of T(E) 
is shown in the Fig. 6. Two peaks in the transmission are associated with the phonon-induced electronic states. 
Their positions and heights depend on the temperature and on the coupling strengths Jo, and w. The important 
feature in the electron transmission is the absence of the peak associated with the resonance state between the grains 
(the bridge site) itself. This happens due to the strong suppression of the latter by the effects of the environment. 
Technically, this peak is damped for it is located at E = where the imaginary part of P{E) reachs its maximum 
in magnitude. To provide the damping of the original resonance the contribution from the environment (including 
the side chain attached to the bridge) to the Green's function (|26|) must exceed the terms Ai,2 describing the effect 
of the grains. This occurs when the inequality 



Ai, 2 < w 2 /^J Q kT (39) 

is satisfied. When the coupling of the bridge to the attached side site is weak, the influence of the environment slackens 
and the original peak associated with the bridge at E = E\ may emerge. At the same time the features originating 
from the phonon-induced states become small compared to this peak. 

So, the effects of the phonons may lead to the damping of the original resonance state for the electron tunneling 
between the metallic islands in the polymer fiber. Instead, two phonon-nduced states appear to serve as intermediate 
states for the electron transport J44 
conductance in molecules Refs. [l 



Environmental induced electron states were discussed in the theory of electron 
4^,13]. For instance, it was shown that low biased current-voltage characteristics 
for molecular junctions with DNA linkers may be noticeably changed due to the occurrence of the phonon-induced 
electron states similar to these discussed in the present Section. 

As we discussed before, one may treat a nanofiber as a set of parallel working channels, any single channel being 
a sequence of grains connected with the resonance polymeric chains. Accordingly, the voltage V applied across 
the whole fiber is distributed among sequental pairs of metallic-like islands included in a single channel. So, in 
practical nanofibers the voltage AV applied across two adjacent grains appears to be much smaller than V (AV/V ~ 
10 _1 — 10~ 3 ). Experiments on the electrical characterization of the polymer fibers are usually carried out at moderately 
high temperatures (T^300K) , so it seems likely that kT > ftw c . Assuming that w ~ lOOmeV, and Jo ~ 204-50meU 
we estimate the separation between the phonon-induced peaks in the electron transmission as 120 -j- 170meV (see 
Fig. 6). This estimate is close to eAU when V takes on values up to 2 4 3 volts. So, the phonon-induced peaks in 
the electron transmission determine the shape of the current- voltage curves even at reasonably high values of the bias 
voltage applied across the fiber. The resulting current-voltage characteristics are shown in the Fig. 7 (left and middle 
panels) . The I — V curves exhibit a nonlinear shape even at room temperature despite the fact that the original 
state for the resonance tunneling is completely suppressed. This occurs because the intergrain transport is supported 
by new phonon-induced electron states. 

It is worthwhile to discuss the temperature dependence of the peak value of the electron transmission which follows 
from the present results. Using the expression ([26)) for the Green's function and the expression (|36[) for P(E) we 
may show that at low temperatures the transmission accepts small values, and exhibits rather weak temperature 
dependence. At higher temperatures (T ~ 100-ftf) the transmission increases fast as the temperature rises and then 
it reduces as the temperature further increases. The peak in the transmission is associated with the most favorable 
conditions for the environment induced states to exist when all remaining parameters (such as Jo and w) are fixed. 
At high temperatures the peaks associated with the environment induced states are washed out, as usual. 

We may compare this result with the temperature dependence of the electron transmission function occurring when 
the bridge between two adjacent grains is directly coupled to the phonon bath. In this case the transmission peak 
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FIG. 6: (Color online) Left panels: Energy dependence of the real (top panel) and imaginary (bottom panel) parts of P(E). 
The curves are plotted at Jo = 20meV, kT S> hoj c . Right panels: The renormalized electron transmission function vs energy. 
The constant Jo equals 20meV (top panel) and 50meV (bottom panel). All curves are plotted assuming w — WOmeV, E\ = 
Ei = 0, T = WOK (dashed lines) and T = 300K (solid lines). 



value may be presented in the form determined by Eqs. (fTTj) . ([24| : 

T(E) = 1 - — ^ (40) 

2(1 + v/l + p 2 ) 2 

The temperature dependencies are shown in the right panel of the Fig. 7. Both curves are plotted at the same value of 
the electron-phonon coupling strength Jo. Comparing them we conclude that at higher temperatures the dependencies 
significantly differ. While the temperature rises, we observe a peak in the electron transmission assuming the indirect 
coupling of the bridge to the phonon bath, and we see the transmission to monotonically decrease when we consider 
the bridge directly coupled to the bath. Correspondingly, we may expect qualitative diversities in the temperature 
dependencies of the current, as well. These diversities originate from the difference in the effects of environment on 
the intergrain electron transport in the cases of direct and indirect coupling of the bridge site to the phonon bath. 
When the bridge is directly coupled to the bath, the stochastic motions in the environment only cause washing out of 
the peak in the electron transmission, and the higher is the temperature the less distinguishable is the peak. However, 
when the bridge is screened from the direct coupling with the phonons due to the presence of the nearby sites, the 
stochastic nuclear motions in the medium between the grains (especially those in the resonance chain) may take a 
very different part in the electron transport in conducting polymers at moderately low and room temperatures. Due 
to their influence, the original intermediate state for the resonance tunneling may be completely suppresed but new 
phonon-induced states may appear to support the electron transport between the metallic-like domains in conducting 
polymer nanofibers. 



V. Conclusion 



Studies of the electron transport in conducting polymers are not completed so far. Several mechanisms are known 
to control the charge transport in these highly disordered and inhomogeneous materials, and their relative significance 
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FIG. 7: The current-voltage characteristics (nA-V) plotted for n — 10, Ai = A2 = 0.5meV, w = WOmeV, Jo = 20meV 
(left panel), and Jo = 50meV (midle panel) at T = 100 A" (dash lines) and T = 300 A" (solid lines). Right panel: Arrhenius 
plot of the peak value of the electron transmission function for Jo = 2077161/, w = WOmeV, Ai = A2 = 0.5meV, and 
To = 100A". The dashed line is plotted assuming the indirect coupling of the bridge to the phonon bath. The solid line is 
plotted assuming that the bridge is directly coupled to the phonons. 



could vary depending on both specific intrinsic characteristics of a particular material (such as crystallization rate 
and electron-electron and electron-phonon coupling strengths) and external factors such as temperature. Various 
conduction mechanisms give rise to various temperature dependencies of the electric current and conductance, which 
could be observed in polymer nanofibers/nanotubes being in conducting state. In the present work we aimed at 
finding out the character of temperature dependencies of both current and conductance provided by specific transport 
mechanism, namely, resonance tunneling of electrons. 

Accordingly, we treated a conducting polymer as a kind of granular metal, and we assumed that the intergrain 
conduction occurred due to the electron tunneling between the metallic-like grains through the intermediate state. 
This scenario for the intergrain electron transport strongly resembles electron transport through molecules/quantum 
dots attached to the conducting leads. In the considered case the metallic-like islands work as the leads and the 
intermediate state in between acts as a single level quantum dot /molecular bridge. Basing on this similarity we did 
apply the well-known Landauer formula to compute the intergrain electron current. This bringed results which agreed 
with experiments on electrical characterization of doped polyaniline-polyethylene oxide nanofibers. 

There are solid grounds to expect significant dissipative effects in the intergrain transport at moderately high 
temperatures. To take into account the effect of temperature we did represent the thermal environment stochastic 
nuclear motions as a phonon bath, and we introduced the coupling of the intermediate site to the thermal phonons. 
As shown in the previous studies of the electron transport through molecules, various dissipative effects may occur 
depending on characteristic features in the interaction of a propagating electron with the environment. Among 
these features we singled out the character of the electron coupling to the dissipative reservoir (phonon bath) as a 
very significant factor. It is likely that in practical conducting polymers both direct and indirect coupling of the 
intermediate state (the bridge) to the environment may occur. 

We did analyze temperature dependencies of transport characteristics for both scenarios and we showed that these 
dependencies differ. Also, wc showed that in general, resonance electron tunneling between the grains results in 
the temperature dependencies of transport characteristics, which differ from those obtained for other conduction 
mechanisms such as phonon-assisted hopping between localized states. Being observed in experiments on realistic 
polymer nanofibers, the predicted dependencies would give grounds to suggest the electron tunneling to predominate 
in the intergrain electron transport in these particular nanofibers. We believe the present studies to contribute to 
better understanding of electron transport mechanisms in conducting polymers. 
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